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(~) I Abstract. We considered the modulational instability of continuous- wave back- 

grounds, and the related generation and evolution of deterministic rogue waves in 
the recently introduced parity-time (7^7')-symmetric system of linearly-coupled non- 

^ , linear Schrodinger equations, which describes a Kerr-nonlinear optical coupler with 

0^ ' mutually balanced gain and loss in its cores. Besides the linear coupling, the over- 

lapping cores are coupled through cross-phase-modulation term too. While the rogue 






[■^^ , waves, built according to the pattern of the Peregrine soliton, are (quite naturally) 



^ 



unstable, we demonstrate that the focusing cross-phase-modulation interaction results 
f^ . in their partial stabilization. For P7~-symmetric and antisymmetric bright solitons, 

CO I the stability region is found too, in an exact analytical form, and verified by means of 

direct simulations. 
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1. Introduction 

It is a generally recognized fact that, independently of the underlying physics, an 
instability of the background is a prerequisite for the emergence of regular or random 
rogue waves (see, e.g., the discussion in Ref. P). In its turn, the instability is 
determined, on the one hand, by the interplay between the dispersion and nonlinearity, 
and, on the other hand, by the competition between losses and gain, if an open system 
is considered. In this latter case, one can speak about dissipative rogue waves [2J, which 
are identified by an enhanced probability of generating high-amplitude pulses. 

In addition to the above-mentioned generic situations, there exist special dissipative 
systems obeying the so-called parity-time (VT) symmetry, i.e., featuring spatially 
separated and exactly balanced gain and loss. These systems are described by non- 
Hermitian Hamiltonians, which may have purely real spectra of eigenvalues, provided 
that the strength of the anti-Hermitian part of the Hamiltonian (which accounts for the 
balanced gain and loss) does not exceed a certain critical value pi H] . 

Optics represents a unifying framework for a variety of wave phenomena. In 
particular, the PT-symmetry was experimentally implemented in coupled optical 
waveguides [5] . Moreover principles of its implementation in plasmonic waveguides [6] 
and in a gaseous mixtures of resonant atoms |7] , were recently proposed. On the other 
hand, optical rogue waves have also been observed in some settings |Hl E] and predicted 
in others, such as periodic arrays of waveguides [10]. While the original ideas of the 
use of the VT symmetry in quantum mechanics imply complex potentials obeying 
condition V{x) = V{—x) [1] (hereafter the overbar stands for complex conjugation), 
in the experimental realization [5] and numerous theoretical studies nonlinear dual- 
core waveguides (couplers), with one core carrying the gain and the other one being 
lossy, were explored as an optical implementation of the "PT-symmetric systems. The 
dual-core systems are described by systems of coupled nonlinear Schrodinger equations 
(NLSEs), one with the gain and the other — with loss. These models and their 
generalizations in a form of sequence of couplers give rise to bright |TTl [121 UHl [E] 
and dark [16] solitons, vortices [H], breathers [17], and describe a switch for solitons 
between the cores [18] . 

As concerns optical rogue waves, there are two major directions of the work in 
this field. The first relates rogue events to the well-known process of super continuum 
generation [191 ES EH [221 123] in optical fibers. The soliton dynamics affects the 
supercontinuum generation process at a very early stage, viz., the fission of higher- 
order solitons [241 [20l [25] . which is followed by multiple interactions of solitons with 
dispersive waves at advanced stages [22l [26l [27]. In particular, the strongest-Raman- 
shifted solitons [28] were proposed as possible candidates for rogue waves [8]. Crests of 
soliton collisions were proposed too, as alternative candidates [29l|30]. Recently, "long- 
lasting" accelerating optical rogue waves with an oblong shape, resembling the shape 
of their oceanic counterparts, were reported [311 [32] . Another approach [HH [HH [35] is 
based on solutions for Akhmediev breathers [36], and, in particular, on the single-peak 
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solution often referred to as the Peregrine soliton (or Peregrine rogue wave) [37], which 
represents a deterministic rogue wave [38] generated by the NLSE recently observed 
experimentally [9]. These works reveal waves which arise from the modulational 
instability (MI) and subsequently disappear, which is consistent with the behavior of 
the famous ship killers in the ocean |39j . 

The main objective of the present paper is to study rogue waves in PT-symmetric 
optical models based on the dual-core couplers. One of our goals is to introduce an 
analog of the Peregrine soliton in this setting. More specifically, we are interested in 
how the presence of the balanced dissipation and gain, i.e., the VT symmetry, affects 
the MI of the background and possibility of the creation of waves localized in space and 
time in such systems. In this context, it is relevant to mention a number of previous 
studies of the deterministic rogue waves carried out in the framework of the coupled 
NLSEs describing two-component matter waves in Bose-Einstein condensates j40j, multi- 
parametric vector solitons, and, in particular, bright-dark- rogue waves jH]. 

The rest of the paper is organized as follows. The model is introduced in Section 
II, which is followed by the analysis of the MI of the continuous-wave (CW) solutions in 
Section III, and the study of rogue-wave solutions, following the pattern of the Peregrine 
soliton, in Section IV. Exact analytical results, verified by direct simulations, for the 
stability of PT-symmetric and antisymmetric solitons in the same system are reported 
in Section V, and the paper is concluded by Section VI. 

2. The model 

We consider a system of linearly coupled NLSEs for field variables ipi and ijj2'- 
. dipi d'^ip 



dz dx"^ 

. ^^l)2 d'^i''. 



r + (xilV^iT + xlH'') ^i + iii^i - i^2, (i) 

+ {x\H^ + xi\H^) i^2 - ^7i^2 - i^i- (2) 



dz dx"^ 

which describes a set of two parallel planar waveguides, with z and x being 
dimensionless propagation and transverse coordinates. Accordingly, the initial-value 
problem corresponds to an optical beam shone into the waveguides input at given z = Zi. 
Alternatively, the model describes a dual-core fiber coupler, where x plays the role of 
the temporal variable [HI [121 US US]. Equations (ll|) and (J2]) are coupled nonlinearly by 
the cross-phase modulation (XPM) ~ x, and linearly by the last terms with respective 
coupling constant scaled to be 1. Lastly, constant 7 > describes the PT-balanced 
gain in Eq. ([T]) and dissipation in Eq. ([2]). In optics, this setting can be realized using 
a system of two lossy parallel-coupled waveguides, doped by gain-providing atoms, in 
which only one waveguide is pumped by the external source of light providing the gain. 
Although the first core carries the gain, its linear coupling to the lossy mate makes 
the zero state in the system neutrally stable, allowing for propagation of linear waves. 
This is the well-known situation, which takes place if the gain/loss term is small enough, 
compared to the linear coupling through which the energy is transferred from the core 
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with the gain to the lossy one, or, more specifically, when 7 < 1 [12] • In such a situation, 
modes can be excited in the system by input beams but do not arise spontaneously. 
Below, without the loss of generality, we restrict the consideration to this case, and 
therefore introduce a convenient parametrization, 

7 = sin(5, 0<5<7i/2. (3) 

Following Ref. pjj, we look for PT-symmetric (+) and antisymmetric (— ) solutions 
to Eqs. ([ID and (^ as 

^2{x,z)=±e^''ilJi{x,z), (4) 

with function ipi obeying the single equation. 



.(9^1 d ipi 



+ {Xi + x)\HS^iT{cos6)tPi. (5) 



dz dx^ 

An observation particularly relevant to the solutions having the form of Eq. (j4]) 
is that the dissipation and gain break the conventional symmetry of the coupler 
The conventional symmetry is now substituted by the following reduction: if 
{ipi{x, z),ip2{x, z)) is a solution of Eqs. ([1]) and ([2]), then pair [ip 2ix, —z),ipi{x, —z)) 
is a solution too. This reduction corresponds to the change 6 ^ n — 6. Therefore, 
below we consider the domain of the variation of 6 to be [0, vr], where values 6 and n — 6 
corresponds to the two different solutions at the same dissipation and gain. In other 
words, intervals < S < 7r/2 and 7r/2 < S < n correspond to the PT-symmetric and 
PT-antisymmetric solutions. 

3. Modulational instability 

Up to a trivial phase shift, CW solutions to Eqs. ([1]) and Eq. ([2]) are 

^j™) = pexp [ikx - ibz + i{-iy5/2] , (6) 

where k represents a background current, and h = k"^ + P^iXi + x) ~ cos 5 (see, e.g., 
[16] for more details), i.e., the amplitudes of the fields are equal in both cores, which is 
natural in view of the necessity to ensure the balance between gain and loss. To study 
the MI of the CW states, we use the standard ansatz, 

Akx—ibz 



^j=P 



^r{-iyS/2 ^ ^^^-^{Pz-n.) ^ ^^^^{Pz-..) 



e 



j = 1,2, with \rij\ lujl ^ 1. Then, two branches f3 = /3i,2(^) of the dispersion relation 
for the stability eigenvalues are given by 



P,{K)^2kK±^ ^K^ + 2p^{Xi + x), (7) 

/32(/t) = 2kK ± ^/[K^ + 2 cos S] [k"^ + 2 cos 5 + 2p2 (^^ _ ^^)] (8) 

We aim to identify parametric domains where the background is subject to the MI. 
Due to the Galilean invariance of underlying Eqs. ([1]) and Eq. (|2]), the instability is 
not affected by boost k. Next, we observe from Eqs. (j?]) and (IE]) that there are three 
different sources of the MI. Firstly, the instability occurs at 

Xi + X<0 (9) 
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This is the "standard" (i.e., observed also for the conservative system of nonhnearly 
coupled NLSEs, without linear coupling) instability stemming from Eq. ([7]) due to 
the long-wavelengths excitations; this domain of the parameters is not influenced by 
gain/dissipation. 

Another instability domain, 

cos5 < max{0,p^(x-Xi)}, (10) 

ensues from Eq. ([8]), and linear coupling between NLSEs gives rise to the appearance 
of this instability domain. Nevertheless, here presence of gain/dissipation [6 ^ 0, vr) 
makes the situation significantly different from that in conservative system (5 = or 
6 = it) j43], as distinct from the previous case. The largest instability growth rate, 
h> = maxK{Im/3(K)}, is 

z/ = Im{/3(K„)} = p'^lxi -xl, K,lr = p'^ix-Xi) -2cos(^, 
in the case 

2cosS<p\x-Xi), (11) 

and 



u 



2y/\[cos5 + p'^iXi - x)]cos5\, ft^ = 



at 



pHx-Xi) <2cos(5<0 and 0<p'^{x-Xi) < 2cos(5 < 2p2(^-^^), (12) 

cf. Eq. (fTTT) . Note that domain (TT2|) disappears in the case of the equal SPM and XPM 
nonlinearities, Xi = X (i-6-, in the PT-symmetric version of the Manakov's system [44j). 
The first consequence of Eq. flTOj) is that for 7r/2 < 5 < vr (antisymmetric solutions) 
the background is unstable irrespective of values of other parameters. The MI regions 
for symmetric solution (0 < 5 < 7r/2) are displayed in detail Fig. [H where the cases of 
focusing {x = ~1) cind defocusing (x = 1) XPM are considered separately. The former 
case [Fig. [11(a)] is the simplest one: here, beyond the fulfillment of condition ([9]) [shown 
by the shadowed region in Fig.[^a)], i.e., at Xi > ~X^ condition fITOl) results in cos 5 < 0, 
i.e., it does not introduce any new domain of the MI. The situation is more complicated 
in the case of the defocusing XPM [Fig.[T]^b)], where along, with xi < —X [the shadowed 
region], there exists another MI domain, generated by Eq. (fTOll . As a result, the CWs 
with large amplitudes, p^ > cos{S)/{2x), are unstable for Xi < X ~ cos((5)/p^. At 
the same time, at 5 — )■ 7r/2, this instability domain approaches the whole area of 

-X < Xi < X- 

Different origins of the MI should naturally lead to different scenarios of its 
development, which we studied by means of direct numerical simulations of Eqs. ([I]), 
([2]). The simulations were performed subject to periodic boundary conditions, and 
with initial excitation of the CW state ([6]) by adding random noise with the amplitude 
amounting to 1% of that of the unperturbed background. 

Starting with the case of weak gain and loss, defined by the symmetric solution 
{6 < 7r/2), in Fig. [2] we show typical results of these simulations for the focusing XPM [see 
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Fig.[2](a), with parameters corresponding to point A in Fig. [T](a)], and for the defocusing 
XPM [see Figs. [2](b)-(d), with parameters corresponding to points B-D in Fig. Wljo), 
respectively]. In panel (a) we observe a "standard" scenario of the development of MI. 
This behavior, being seemingly expectable, nevertheless reveals a noteworthy feature of 
the PT-symmetric system, which behaves similarly to its Hamiltonian counterpart (at 
least, in significant initial intervals of the propagation). In particular, we observe that 
the power is distributed between the two waveguides. 

The situation changes significantly when one consider the defocusing XPM, even 
if Eq. (Q is satisfied, i.e., the MI has the same nature as in the conservative system. 
Indeed, in Fig. [2](b) we observe a rather fast power transfer from the lossy waveguide 
to the one with the gain, accompanied by fast growing peaks. Obviously, such peaks 
can be described by a single NLS equation (^ with ip2 = 0. The observed behavior is 
due to the focusing SPM, xi, and therefore is not significantly altered even when one 
passes from the domain of parameters ([9]) [Fig. [2](b)] to the one defined by Eq. (fTOl) . as 
shown in Fig. M^c). A significant change, i.e., the third scenario of the evolution of the 
MI, appears when the SPM is defocusing too [Fig. [2](d)]. This is the case where the MI 
occurs only due to the imbalance of the gain and loss, resulting in nearly homogeneous 
grow (decay) of the field in the waveguide with gain (dissipation), respectively. 

Examples of the modulational instability and stability for the CW solution with 
nonzero wave vector k (current) are presented in Fig. [3l Here we restrict our 
consideration of the MI with the focusing SPM, xi < 0, but when Xi + X < [point C 
in Fig. [U^b)]. Thus, the evolution of the MI in this case occurs according to the same 
scenario as for k = 0, cf. Figs M^c) and Figl3]^a). At the same time, the respective MI 
peak is shifted in the positive direction of the x-axis, which coincides with the direction 
of the current. Meanwhile, in the domain where the CW state is predicted to be stable 
[above the green line, in Fig{Tt^b) — e.g., at point C ], the stability is confirmed by the 
numerical simulations, see Figj3]^b). 
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Xi=-X 


(a) 
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1 2 4 


P 
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Figure 1. Domains of tlie MI: sliaded regions in both panels, as well as the region 
under the respective curve, corresponding to different gain/loss coefficient (S), in panel 
(b), in the (xi, /o)-plane for different 5 (as indicated in the panels) and for fixed XPM 
coefficients, x = — 1 (a)i or x = 1 (b)- Capital letters indicate parameters chosen for 
displaying the evolution in figures following below. 
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Figure 2. The evolution of field components \ipi{x, z)p and |V'2(2;, z)\'^ (left and right 
columns) of the plane-wave solution with parameters k = 0, S = 7r/4, p — 1.604, 
XI = 0.5, X = -1 (a), P = 0.76, xi = -1-5, X = 1 (b), P = 0.79, xi = -0.5, x = 1 
(c) and p — 0.98, xi — 0.25, x — ^ (d). Parameters of panels (a), (b), (c) and (d) 
correspond to points A, B, C and D in Fig. [1] respectively. 



4. The Peregrine soliton in PT-symmetric system: the case of xi + x < 

Turning now towards studying the Peregrine soliton (rogue wave) propagating against 
an unstable background we start with the case ([9]). This readily allows one to write 
down the Peregrine solution of Eqs. (Ill), (E]) in the form (j = 1, 2) ^\ H5] 

^,-(x,z)=pe("^)''^/2+'^"-'''^x 

_4 1 - 2t (xi + x) P^z 1 

l-2{xi + x)p'{x-2kzY + A{xi + xfp^zA' 
Notice that, when \z\ — )■ oo, or, equivalently, |x| — )■ oo, solution (IT3l) merges into the 
background given by Eq. ([6]). Below, we separately consider two cases: the Peregrine 
soliton, based on the background without the current {k = 0), and current-based 
Peregrine solution, with k ^ 0. 

Examples of the Peregrine solutions whose backgrounds (without the current. 
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Figure 3. The evolution of field components |^i(a;,z)p |^2(a;, 2)^ (left and right 
columns) of the plane-wave solution with parameters k — 0.2, S ~ 7r/4, xi — —0.5, 
X = 1, p = 0.79 (a), or p ~ 0.4 (b). Parameters of panels (a),(b) correspond to points 
C, C in Fig. [TJb), respectively. 



k = 0) correspond to points A and B in Fig. [H are depicted in Fig. HI The spatial 
evolution of \il!j{x,z)\'^ was obtained by the numerical simulations of Eqs.(IT]), (|2]) with 
the initial condition corresponding to the Peregrine soliton ( IT3l) at z = Zi = —4 [Figs. 
m^a) andlH^b)], or z = Zi = —2 [Figs. Ht^c) andHfd)]. In the case of the defocusing 
SPM and focusing XPM [FigsJU^a) andlD^b)], the central peak, corresponding to the 
Peregrine solution, appears at x = 2; = 0, before MI peaks. At the same time, in the 
case of the focusing SPM and defocusing XPM [Figs. Hl^c) andlD^d)], the appearance 
of the Peregrine-soliton peak a.t x = z = causes further growth of the peak in the 
first (gain-pumped) component, and decrease in the second (lossy) core. Notice that 
the structure of the rogue-wave evolution in this case resembles the respective scenario 
of the MI development for the same parameters, as suggested by the comparison of 
Figs. El^b) and Figs.lD^c). 

In the case when the background carries the current {k ^ 0), the central peak of 
the Peregrine solution moves with group velocity 2k in the positive direction of x-axis, 
as seen in Eq.( TT3|) and confirmed by Figs. |5]^a,b). Also for the focusing-XPM {x < 0) 
case, the PT-symmetric {5 < 7r/2) rogue wave is more "stable" (in the sense that the 
MI peaks appear after at a longer propagation distance after the principal rogue-wave 
peak), see FigJSt^a), if compared to the PT-antisymmetric wave with 6 > tt/2, see Fig. 

Eh). 

In order to describe this rogue wave "stability" quantitatively, we will use one of 
the principal properties of Peregrine solution, which follows from Eq. (fT3|) . namely 



V^, 






ips-jixjz). If the phase of the solution is not taken into account, this 



property turns into \iIjj{—x 



S 



00 



|^j(a;, z)p. Thus, we introduce the discrepancy as 



21 2 



i{-x,-Zi)\ +\ij2{-x,-Zi)\ -\tlJi{x,Zi)\ -\tlj2{x,Zi)\] dx 
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\W 



4y\,/^a) J|^/ 4_/X^^^a) 




Figure 4. Peregrine solutions in tlie T'T-limit for p — 1.604, xi = 0-5, X = ^li 
S = 7r/4 (a), or S ^ Sn/A (b); p ^ 0.76, xi = -1-5, x == 1, -5 == 7r/4 (c), or (5 = 37r/4 
(d) . Parameters of panels (a) and (b) correspond to the point A, while those of panels 
(c) and (d) - to point B in Fig[T] 



in order to eliminate phase effects. In the ideal case, where the shape of the rogue 
wave coincides with the Peregrine solution flT3|) . the discrepancy is zero, S = 0. Thus, S 
serves to measure how much the numerically obtained solution differs from the Peregrine 
soliton, or in other words, how much the chaotic nature of MI influences the Peregrine 
solution. The results are depicted in Figs. Et^c) and Et^d). For the focusing XPM 
[Figj5]^c)] and for 6 < 7r/2 the discrepancy abruptly grows at p > 1. In the same time, 
in this range of 6 the discrepancy almost does not depend on 6 [the lines for 6 = tt/A, 
6 = 97r/20, and 6 = ll7r/20 are indistinguishable on the scale of Fig. EJ^c)]. Meanwhile, 
for 6 > Tx 12 the situation is opposite: the discrepancy increases with b (compare the 
lines for b = ll7r/20 and 6 = Stt/A). For the defocusing XPM [Fig. Et^d)], discrepancy 
S decreases with the increase of 6 in the whole range of < 5 < tt, while an abrupt 
growth happens at p > 0.4. As a result, for the focusing XPM, the PT-symmetric rogue 
wave with 6 < it/2 is more stable than the its antisymmetric counterpart, while for the 
defocusing XPM the situation is opposite. 
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""0 0.3 



Figure 5. (a,b) Current-based Peregrine solutions in the PT-symmetric system for 
k = 0.6, p ^ 1.0, xi = 0.5, X = -1, and S = 7r/4 (a), or S = Sn/A (b); (c,d) 
Discrepancy S vs amplitude p for different gain/loss coefficients 5 (indicated in panels) 
and for xi = 0-5 and the focusing XPM, with x = ^1 (c) or xi = —1.5 and the 
defocusing XPM, with x — ^■ 



5. Bright solitons 

Obvious bright-soliton solutions of Eq. ([5]) with arbitrary amplitude tj are available too, 
for xi + X < 0: 



^j 



V 



exp 



-iy- + cos6 + -r]^ 



(14) 



A/|Xi + x|cosh {r]x/V2) 
j = 1, 2, where, as above, intervals < S < 11/2 and 7r/2 < 6 < tt correspond for the 



PT-synimetric and antisymmetric solitons, respectively. Using results from Refs. [TT] 
and [47J, an exact stability boundary for the symmetric and antisymmetric solitons, 
against small perturbations breaking the respective symmetry or antisymmetry, can be 
predicted in the following analytical form: 

16 (-Xi -x)cos5 



vl 



(15) 



(V-25xi + 7x - 3V-X1 - X) (V-25xi + 7x + V'Xi - x) 
the solitons being stable at 77^ < rj"^^. 

This result makes sense when Eq. flTSj) yields a positive value, otherwise the 
PT-symmetry-breaking bifurcation does not occur, and the stability may only be 
studied numerically [in addition to the instability mode represented by Eq. ( !T5|) . other 
instabilities are possible too]. In particular, condition (ITSll cannot simultaneously hold 
for the PT-symmetric and antisymmetric solitons. Further, because the existence of 
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Figure 6. (a) Stability boundaries for T'T-symmetric and antisymmetric solitons (|T4)) 
are shown in the plane of the gain-loss coefficient 7 (recall 7 = sin J) and soliton's 
amplitude 77, in the case of opposite signs of the SPM and XPM coefficients. The 
boundaries for the PT-symmetric (SI) and antisymmetric (Al) solitons are shown by 
solid blue and red lines, respectively, for x = 1 and xi — —1-5. The dotted black curve 
labeled SI (an) displays the analytical counterpart of the SI boundary, as predicted by 
Eq. (jisp . The stability boundary for antisymmetric solitons, depicted by the dashed 
green curve (A2), pertains to x = ^1 and xi = 0.5. (b) The stability boundaries 
in case of identical signs of the SPM and XPM coefficients, x = ~1 and xi = ^3. 
The boundaries (S3) and (A3) for symmetric and antisymmetric solitons are shown by 
solid blue (S) and dashed red curves (A), respectively. The dotted black curve labeled 
S3(an) is the analytical counterpart of the latter boundary, as predicted by Eq. ([T5|) . 



the solitons of either type requires xi + X < 0, the condition of 77^ > actually holds 
for the PT-synimetric solitons at — Xi > —Xj ^^^ for the PT-antisymmetric solitons — 
in the opposite case, at —Xi < — X- 

We have performed direct simulations of the evolution of perturbed solitons within 
the framework of Eqs. ([1]) and ([2]), aiming to identify stability borders for the VT- 
symmetric and anti-symmetric solitons, and, in particular, to verify the analytical 
prediction ( [T5|) . Perturbations were introduced by adding 2% to the amplitude of 
component, and reducing 2% from the other. Figures |6]^a) and Mjo) display the so 
identified stability boundaries in the cases of opposite and identical signs of x and xi, 
respectively. For the sake of comparison with Ref. [HI [12], we demonstrate these borders 
as a function of 7 , rather than S [see Eq. ([3])]. 

The numerically found stability boundaries are close to their analytical 
counterparts. Some discrepancy between them is explained by the fact that some 
solitons, which are stable against infinitesimal perturbations, may be destabilized by 
finite-amplitude excitations. 

Typical examples of the unstable and stable evolution of antisymmetric solitons, 
taken on both sides of the stability boundary, are demonstrated in FigJTJ for x = —1 
and Xi = 0-5. The quick stabilization of the symmetric soliton in the same case is 
demonstrated in Figl8]for a large amplitude, rj = 3. Actually, the PT-symmetric solitons 
are stable for all rj in this case, the stability border being relevant for the antisymmetric 



ones. 
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Figure 7. Stable and unstable evolution of both components of antisymmetric solitons 
([H| at X = —1 J Xi = 0-5 and S — tt — arcsin(0.2). Panels (a) and (b) pertain to the 
unstable dynamics of the soliton with amplitude i] = 0.6, while panels (c) and (d) 
pertain to the stable soliton with rj = 0.3. 




Figure 8. Quick stabilization of a symmetric soliton with a large value, r/ — 3, for 
X — —1, Xi — 0.5 and 7 = 0.2. Dynamics of the u-component is similar to that shown 
here for the w-component. 



It is relevant to note too that, in the Manakov's hmit, Xi = X P^? the stabihty 
boundary predicted by Eq. f lTSj) diverges. Indeed, direct simulations demonstrate that 
all the solitons are stable in this case. 



6. Conclusion 

In this work we have considered the MI (modulational instability) of CW backgrounds 
and the emergence and evolution of rogue waves in the system of linearly-coupled 
"PT-symmetric coupled NLSEs. We have shown that the focusing XPM nonlinear 
interactions extend the effective stability region for the rogue waves of the Peregrine's 
type. The system can support nondissipative rogue waves too. The stability region for 
PT-symmetric and antisymmetric solitons was found in the exact analytical form and 
verified by direct simulations. It may be interesting to extend the analysis for (2D) 
versions of the system, which may have realizations in nonlinear optics, cf. Ref. 
and references therein. 
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